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Abstract 

We present a QCD calculation of the transverse momentum distribution of photon pairs produced 
at hadron colliders, including all-orders soft-gluon resummation valid at next-to-next-to-leading 
logarithmic accuracy. We specify the region of phase space in which the calculation is most reliable, 



compare our results with data from the Fermilab Tevatron, and make predictions for the Large 
Hadron Collider. The uncertainty of predictions for production of diphotons from fragmentation of 
final-state quarks is examined. 
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Introduction. A Higgs boson with mass between 115 and 140 GeV may be identified at 
Preprint submitted to Elsevier Science March 3, 2006 



hadron colliders through its decay into a pair of energetic photons, a challenging prospect at 
the Large Hadron Collider (LHC) in view of the intense background from hadronic produc- 
tion of non-resonant photon pairs [1] . Theoretical predictions of these background processes 
may be of substantial value in aiding search strategies. Moreover, the perturbative quantum 
chromodynamics (QCD) calculation of photon-pair production is of theoretical interest in 
its own right, and data from the Tevatron collider offer an opportunity to compare and test 
results against experiment. 

In this paper, we present a new calculation of the diphoton cross section in perturbative 
QCD. We include contributions from all perturbative subprocesses (quark-antiquark, quark- 
gluon, antiquark-gluon, and gluon-gluon) to next-to-leading (NLO) accuracy. In addition, to 
describe properly the behavior of the transverse momentum Q T distribution of the pairs in 
the region in which Qr < Q, where Q is the invariant mass of the photon pair, we include the 
all-orders resummation of soft and collinear logarithmic contributions up to next-to-next- 
to-leading log (NNLL) accuracy. This calculation goes beyond the previous resummation 
treatments of diphoton production [2,3]. Its components are summarized briefly below, and 
a more complete discussion is presented elsewhere [4]. 

A full treatment of photon pair production requires that we address the contributions from 
non-perturbative processes, such as it and rj meson decays, and the quasi-collinear fragmen- 
tation of quarks and gluons into photons. Elaborate isolation procedures are applied by the 
experiments to reduce these long-distance contributions, procedures that are only approx- 
imately reproducible theoretically. Some final-state fragmentation contributions invariably 
survive the isolation, especially at the LHC, where the efficiency of isolation is reduced by 
event pile-up and the large number of energetic hadronic fragments in each event. A new 
feature of diphoton production, with respect to single photon production, is the prospect 
that both photons may be produced from fragmentation of the same final-state parton. This 
fragmentation contribution is expected to be most influential in the region in which both 
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the diphoton invariant mass and the separation A(p between the azimuthal angles of the two 
photons are relatively small, Q < Qt and A<p < ir/4. 

Diphoton production is characterized by large radiative corrections, distributed in a complex 
pattern over the accessible phase space. The influence of initial-state gluon radiation on the 
predicted transverse momentum distributions can be evaluated to all orders with the Collins- 
Soper-Sterman (CSS) resummation procedure [5], the method that we follow. Our results 
are implemented in a Monte-Carlo integration program RESBOS. We use a simple, efficient 
approximation for the fragmentation contributions. We compare our results with data from 
the Collider Detector at Fermilab (CDF) collaboration at pp collision energy y/S = 1.96 TeV 
and integrated luminosity 207 pb _1 [6], and we observe good agreement. We make several 
suggestions for a further more differential analysis of the data that would allow refined tests 
of our calculation. In view of theoretical uncertainties associated with the fragmentation 
component of the cross section, and the presence of other large radiative corrections, we 
question the conclusion in Ref. [6] that the inclusion of single-photon fragmentation con- 
tributions within the NLO calculation of Ref. [7] uniquely explains the observed kinematic 
distributions of the diphotons at the Tevatron. We also include predictions for diphoton 
production at the LHC. 

Analytical Calculation. The CSS resummation method is used in Refs. [2,3] to treat the 
direct production of photon pairs from qq : ( qg, and gg scattering. The NLO perturbative 
cross sections (i.e., cross sections of 0(a s ) in the qq and qg channels [8,9,10], and 0{a z s ) in 
the gg channel [3,11,12,13]) are included as a part of the resummed cross section. Singular 
logarithms arise in the NLO cross sections when the transverse momentum Q T of the 77 
pair is much smaller than its invariant mass Q. These logarithms are resummed into a 
Sudakov exponent (composed of two anomalous-dimension functions A(fj) and B(/i)) and 
convolutions of the conventional parton densities f a (x,fj,p) with Wilson coefficient functions 
C. In Refs. [2,3], the functions A(/i), B(/i) and C are evaluated up to order a 2 s1 a s , and 
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a s , respectively. An approximate expression is used there for the C-function of order a s in 
the gg subprocess (borrowed from the gg — > Higgs resummed cross section). In this work, 
we include the exact C-function of order a s for gg — > 77X [14] and 0{a 2 s ) expressions for 
A(fi) and B(/jl) in all subprocesses [14,15,16]. These enhancements elevate the accuracy of the 
resummed prediction to the NNLL level. We use an improved model for the non-perturbative 
contributions at large impact parameter [17]. When expanded in a series in a s , the resummed 
predictions for the total rate, 77 invariant mass, and 77 rapidity (y) distributions are equal 
to the fixed-order QCD cross sections, augmented by higher-order contributions from the 
integrated Q T logs. The resummed Q T distribution is well-behaved as Q T — > 0, unlike its 
fixed-order counterpart which is singular in this limit. As Qt grows, our resummed cross 
section crosses the perturbative NLO cross section at Qt ~ Q, and, for each Q and y, we 
switch from the resummed to the NLO cross section for values of Q T above this point. 

A fragmentation singularity arises in the matrix element when the momentum of a photon 
is collinear with that of an outgoing quark or gluon. The fragmentation singularities do not 
appear in the resummed terms since those correspond to initial-state radiation. At the lowest 
order, the fragmentation singularity appears in the qg — > 775 channel and is proportional 
to P^gi^z) I [n — 4) in n-dimensional regularization, where P 7< _ 9 (z) is the q — > 7 splitting 
function, and z is the fraction of the fragmenting quark's light-cone momentum carried by 
the photon. The fragmentation singularity is subtracted from the direct contribution. It is 
resummed in the photon fragmentation function D y (z) through the introduction of a "one- 
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fragmentation" contribution q + g — > (q — ► 7 X) +7, where u (q — ► 7X)" denotes collinear 
production of a photon from a quark. For a wide class of two-to-two partonic processes, 
such as qq — > qq, etc., there is a second type of "one- fragmentation" contribution that 
arises in low- mass photon-pair production (Q < Qt)- In this case, a final-state parton may 
fragment into a low-mass pair of photons, a process described by a different fragmentation 
function D 77 (zi, z-i). This new contribution is not included in the existing calculations. "Two- 
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fragmentation" contributions arise in processes like g + g -> (q ^ 1 X) + (q ^ 7 X) and 
involve convolutions with two functions D y (z) (one per photon). 

Isolation constraints must be imposed on the inclusive photon cross sections before the com- 
parison with data. Isolation can be applied to the cross sections at each order of a s [18,19,20]. 
The magnitude of the fragmentation contribution is controlled by the isolation procedure 
chosen and can be strongly affected by tuning the quasi-experimental isolation model. An 
isolation condition in a typical measurement requires the hadronic activity to be minimal 
(e.g., comparable to the underlying event) in the immediate neighborhood of each candi- 
date photon. Candidate photons may be rejected because of energy deposit nearby in the 
hadronic calorimeter, which introduces dependence on the calorimeter cell geometry, or be- 
cause hadronic tracks are present near the photons. A theory calculation may approximate 
the experimental isolation by requiring the full energy of the hadronic remnants to be less 
than a threshold "isolation energy" E™° in the cone AR = y/ Arj 2 + Ap 2 around each pho- 
ton, with Ar] and Aip being the separations of the hadronic remnant (s) from the photon 
in the plane of pseudo-rapidity i] and azimuthal angle ip. The two photons must also be 
separated in the i] — ip plane by an amount exceeding the approximate angular size Ai? 77 of 
one calorimeter cell. The values of E™°, AR, and A_R 77 serve as crude characteristics of the 
actual measurement. The size of the fragmentation contributions depends tangibly on the 
assumed values of E™°, AR, and AR~ n , as is shown below. 

We find it sufficient in our work to use a simplified fragmentation model to represent the 
isolated cross section. We regularize the fragmentation region by imposing a combination of 
a sharp cutoff E™ ° on the transverse energy E T of the final-state quark or gluon and smooth 
cone isolation [21]. We impose quasi-experimental isolation by rejecting an event if (a) the 
separation Ar = — rj^) 2 + ((/? — (/? 7 ) 2 between the final-state parton and one of the pho- 
tons is less than AR, and (b) Et of the parton is larger than E^°. This condition excludes 
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the singular fragmentation contributions in the finite-order qg cross section at Ar < AR and 
E T > Et°- The fragmentation contributions at Ar < AR and E T < E^° are suppressed by 
rejecting events in the AR cone that satisfy E T < x(Ar), where x(Ar) is a smooth function 
satisfying x(0) = 0, x(AR) = E^°. This "smooth cone isolation" [21] transforms the non- 
integrable fragmentation singularity associated with D 1 {z) into an integrable singularity of a 
magnitude dependent on the functional form of x(Ar). Infrared safety of the cross sections is 
preserved as a result of smoothness of x(Ar). The cross section for direct contributions is ren- 
dered finite by this prescription without the explicit introduction of fragmentation functions 
Dj(z). For our smooth function, we choose x(Ar) = E^°(l — cos Ar) 2 / (1 — cos AR) 2 . Modi- 
fications to the function x(Ar) l e& d to only mild variations of our predicted Qt distribution 
for Q T < E^°. 

In our calculation, we use the electroweak parameters [22] Gf = 1.16639 x 1CT 5 GeV~ 2 , 
m z = 91.1882 GeV, and m w = 80.419 GeV. We use two- loop expressions for the running 
electromagnetic and strong couplings and as(/x), as well as the NLO parton distribution 
function set CTEQ6M [23] and set 1 of the NLO photon fragmentation functions from 
Ref. [24]. Our choices of the renormalization and factorization scales are the same as in 
Ref. [2]; in particular, we set Hr = Hf = Q in the finite-order perturbative expressions. 

In impact parameter (b) space, used in the CSS resummation procedure, we must integrate 
into the non-perturbative region of large b. Contributions from this region are known to be 
suppressed at high energies [25], but it is important nevertheless to evaluate the expected 
residual uncertainties. We use a model for the non-perturbative contributions ("revised b* 
model") based on the analysis of Drell-Yan pair and Z boson production in Ref. [17]. A non- 
perturbative Sudakov function for the factorization constant C 3 = 2e~ 7E 1.123 is used 
here to describe the non-perturbative terms in the leading qq — > 77 channel [17]. We neglect 
possible corrections to the non-perturbative contributions arising from the final-state soft 
radiation in the qg channel, as well as additional v^S* dependence affecting Drell-Yan-like 
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processes at x < 10~ 2 [26], as those exceed the accuracy of the present measurements at 
the Tevatron. The non-perturbative function in the gg — > 77 channel is approximated by 
multiplying the non-perturbative function for the qq channel by the ratio Ca/Cf = 9/4 of 
the color factors Ca = 3 and Cf = 4/3 for the leading soft contributions in the gg and qq 
channels. Comparing our results based on the "revised 6* model" with those obtained with 
the original 6* approach, we find at most 10% differences in our predicted da/dQx at the 
lowest values of Qt at the Tevatron collider energy, and smaller differences at larger values 
of Qt, all well within the experimental uncertainties. The differences are even smaller at the 
LHC energy [4] . 

Comparison with Tevatron Data. Our analysis provides a calculation of the triple-differential 
cross section da / ' dQdQTdAip . Its relevance is especially pertinent for the transverse momen- 
tum Qt distribution in the region Qt < Q, for fixed values of diphoton mass Q. It would be 
best to compare our multi-differential distribution with experiment, but the published col- 
lider data tend to be presented in the form of single-differential distributions in Q, Qt, and 
Aip, after integration over the other variables. We follow suit in order to make comparisons 
with the Tevatron collider data, but we comment on the features that can be explored if 
more differential studies are made. In accord with CDF, we impose the cuts |y 7 ] < 0.9 on 
the rapidity of each photon, and p T > p T min = 14 (13) GeV on the transverse momentum 
of the harder (softer) photon in each 77 pair. We choose = 1 GeV, AR = 0.4, and 
Ai? 77 = 0.3, unless stated otherwise. 

The invariant mass (Q) distribution is shown in Fig. la. It exhibits a characteristic lower 
kinematic cutoff at Q ps P^minP Train ~ 27 GeV. Our calculation (RESBOS) agrees well 
with the data. In this figure we also show the perturbative QCD contributions evaluated 
at finite order, represented by the DIPHOX code [7]. Unless specified otherwise, the scales 
= fJ-F = Q are used to obtain the DIPHOX results presented here. The overall agreement 
between the two calculations is anticipated, since both evaluate the inclusive rates at NLO 
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pp -> yyX, CDF Run-2, 207 pb 1 PP ~> YYX, CDF Run-2, 207 pb 
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(a) (b) 

Fig. 1. (a) Invariant mass and (b) transverse momentum distributions of diphotons. Data from 
the CDF Run-2 measurement [6] are compared to our calculation (RESBOS) and the DIPHOX 
calculation. 

accuracy. The differences are due to different isolation prescriptions, resummation of higher- 
order logarithms as well as NLO contributions to the gg channel in our calculation, and 
single-photon one- and two-fragmentation contributions in DIPHOX. 

The transverse momentum (Qt) distribution of diphotons is shown in Fig. lb. The finite- 
order calculation, represented here by DIPHOX, displays an unphysical logarithmic singu- 
larity as Qt — * 0. In our work, the initial-state small-Qr singularities are resummed in the 
CSS formalism, resulting in a reasonable overall shape of the cross section at any Qt- The 
fragmentation contributions exhibit a double-logarithmic singularity when Qt approaches 
E%° from below [7], as it is evident in the DIPHOX Q T distribution for J5j?° = 4 GeV. No 
such singularity is present in our Qt distribution, which instead has a mild discontinuity at 
the point Qt = where we switch from the quasi-experimental to smooth-cone isolation. 
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pp -> yyX, CDF Run-2, 207 pb"' ,_, pp -> yyX, Vs = 1 .96 TeV 
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(a) (b) 

Fig. 2. (a) Distribution over the azimuthal separation Ay? between the two photons. Data from 
the CDF Run-2 measurement [6] are compared with our calculation (RESBOS) and the DIPHOX 
calculation for different isolation parameters; (b) same as (a), with an additional cut Qt < Q on 
the diphoton momentum. Our cross sections are evaluated with the factorization scales [If = Q 
(lower curve) and 0.5 Q (upper curve) in the finite-order contribution. 

For the same value = 1 GeV, our distributions and those of DIPHOX agree well at 
large Qt, as a result of our smooth matching of the resummed cross section to the NLO 
cross section. In the two highest-Qy bins, the CDF central values lie above the two theory 
predictions. While the observed excess of events in this "shoulder" region is not significant 
compared to the present experimental errors, it has been discussed as a possible indication 
of enhanced fragmentation contributions in the Tevatron data [6,27]. 

The parameters in DIPHOX can be adjusted to bring its results into agreement with the 
data in the shoulder region (cf. the dash-dot curves in Figs, lb and 2a). The cross section 
in that region is enhanced if a smaller factorization scale is used, and if the isolation energy 
Elf is increased. The dash-dot curves in Figs, lb and 2a are obtained with /ip = 0.5 Q and 
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Et° = 4 GeV, compared to the nominal value of E™° = 1 GeV in the CDF publication. 
In the shoulder region, the increase in to 4 GeV strongly enhances the DIPHOX cross 
section to the value shown in the CDF publication. The magnitude of the one-fragmentation 
cross section associated with D y (z) is increased on average by 400% when E™° is increased 
from 1 to 4 GeV. 

Our calculations show that most of the shoulder events populate a limited volume of phase 
space characterized by Acp < 1 rad, Q < 27 GeV, and Qt > 25 GeV. The location of the 
shoulder in the Q T distribution is sensitive to the value of the cut on the minimum transverse 
momentum, p T , of the individual photons, moving to larger Qt if these cuts are raised. It has 
also been noted [27] that non-zero values of and p^ disallow contributions with small Qt 
if the azimuthal angle separation between the two photons is small, Aip < ir/2. The excess 
of the experimental rate over our prediction in the region Atp < 0.6 radian (cf. Fig. 2a) 
contributes the bulk of the excess seen in the shoulder in the Q T distribution in Fig. lb. We 
note, in addition, that the excess at small Aip and large Qt is characterized by Qt > Q- 

From a theoretical point of view, when Q T > Q, as in the shoulder region, the calculation 
must be organized in a different way [28,29] in order to resum contributions arising from the 
fragmentation of partons into a pair of photons with small invariant mass. In addition, a small 
azimuthal separation Aip often implies that the photons are produced at polar angles 9* ~ 
or 7r in the Collins-Soper diphoton rest frame [30]. The matrix element for the Born scattering 
process qq — > 77 diverges as |cos0*| — > 1. Large QCD corrections are known to exist when 
I cos 6**| ~ 1 at any order of the strong coupling strength a s . Radiation of additional partons 
at higher orders regularizes the singularity of the quark propagator, yet the enhancement 
of the cross section is still felt at large |cos#*|. At small Q T , the |cos#*| 1 contributions 
are excluded by the cuts p T > 14 (13) GeV on the transverse momenta of the individual 
photons. If, however, the diphoton system is boosted in the transverse direction (Qt > Q), 
contributions with |cos0*| ~ 1 and substantial rapidity separation |y 71 — y 12 \ > 0.3 are 
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allowed in the event sample. 



Adequate treatment of the light 77 pairs and large- 1 cos 6**1 contributions is missing in both 
our calculation and DIPHOX. The presence of higher-order contributions is reflected in 
the sensitivity of the DIPHOX prediction at small Aip to variations of E™°, factorization 
scales, and the angular separation Ai? 77 between the photons [27]. In view of the theoretical 
uncertainties in the calculation of the fragmentation contributions, and the likely presence 
of other types of radiative corrections, we suggest that more theoretical and experimental 
effort is needed to firmly establish the origin of the excess rate in the CDF data at large 
Qt and small Atp, whether from single-photon fragmentation as implemented in DIPHOX 
or/and other types of enhanced scattering contributions. 

The theoretical ambiguities arise in a small part of phase space, where the cross section is 
also small. Our theoretical treatment is most reliable in the region in which Qt < Q- When 
the Qt < Q selection is made, the contributions from Aip < tt/2 are efficiently suppressed, 
and dependence on tunable isolation parameters and factorization scales is reduced (cf. 
Fig. 2b). The fixed-order predictions agree well between our calculation and DIPHOX, while 
our resummmed cross section also provides an accurate description of the rate at small values 
of Qt- After the selection Q T < Q, we expect that the large Qt shoulder will disappear in 
the experimental Qt distribution. 

An important prediction of the resummation formalism is a logarithmic dependence on the 
diphoton invariant mass Q. In Fig. 3a, we show the resummed transverse momentum distri- 
butions for various intervals of Q. The Qt distribution is predicted to broaden with increasing 
Q. The average values of Q T are (Qt) = (6.5, 8.1, 10.7, and 12.6) GeV for invariant masses in 
the intervals (30-35, 35-45, 45-60, and 60-100) GeV, respectively. To compute these averages, 
we integrate over the range Q T = to 200 GeV. We urge the CDF and DO collaborations to 
verify this predicted broadening with Q. 
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pp -> yyX, Vs = /. 96 TeV pp -> yyX, Vs = 14TeV 
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(a) (b) 

Fig. 3. Resummed transverse momentum distributions of photon pairs in various 77 invariant mass 
(Q) bins at (a) the Tevatron and (b) the LHC. The curves are calculated with the cuts specified in 
the text. The cross sections are normalized to the total cross section in each bin of Q. We note that 
our predictions are most reliable in the region Qt < Q, but we plot the curves over the full range 
of Qt, using the procedure described in the text to switch from the resummed to the finite-order 
perturbative results for Qt > Q- 

Results for the LHC. To obtain predictions for pp collisions at yS = 14 TeV, we employ 
the following cuts on the kinematics of the individual photons. For each photon, we require 
transverse momentum p T > 25 GeV and rapidity |y 7 | < 2.5. We impose a somewhat looser 
isolation restriction than for the Tevatron study, requiring less than = 10 GeV of extra 
transverse energy inside a cone with AR = 0.7 around each photon. 

Figure 3b shows the resummed transverse momentum distributions for various selections of 
diphoton invariant mass at the LHC. The plot shows the broadening of the Qt distribution 
with increasing mass: in the ranges (55-65, 65-95, 95-130, and 130-250) GeV the values of 
(Qt) are (14, 17, 25, and 33) GeV. At the LHC, we integrate from Q T = to 250 GeV 
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to obtain the averages. For the mass range appropriate in the search for a Standard Model 
Higgs boson, e.g., 115 to 130 GeV, the diphoton background that we consider in this paper 
has (Qt) ~ 27 GeV, to be compared with the expectation for the signal of ~ 40 GeV [25]. 
The harder transverse momentum distribution for the signal arises because their is more soft 
gluon radiation in the dominant gluon- fusion Higgs boson production process [25]. Additional 
predictions for the LHC are presented in Ref. [4]. 

Summary. We present a new QCD calculation of the transverse momentum distribution of 
photon pairs produced at hadron colliders, including all-orders resummation of initial-state 
soft-gluon radiation valid at next-to-next-to-leading logarithmic accuracy. This calculation is 
most appropriate for values of 77 transverse momentum Qt not in excess of the 77 invariant 
mass Q. Resummation changes both the shape and normalization of the Q T distribution, 
with respect to a finite-order calculation, in the range of values of Qt where the cross 
section is largest. Comparison of our results with data from the Fermilab Tevatron shows 
good agreement, and we offer suggestions for a more differential analysis of the Tevatron 
data. We also include predictions for the Large Hadron Collider. 

Our calculation accounts for the effects of soft gluon radiation on transverse momentum 
distributions through all orders of a s . The NLO calculation with inclusion of single-photon 
fragmentation [7] is another important approach to 77 production. However, theoretical 
uncertainties are present in the rate of fragmentation contributions associated with the kine- 
matic approximations and tunable parameters in the quasi-experimental isolation condition. 
For Qt > Q (Aip < ir/2), new types of higher-order contributions are expected to enhance 
the rate above our predictions. The interpretation of the region of small Aip remains am- 
biguous, as several distinct processes may contribute to the enhanced rate. This interesting 
region warrants further theoretical investigation. With the contributions from the Qt > Q 
region removed, our calculation describes the leading contributions in the qq + qg and gg 
diphoton production channels at NNLL accuracy. 
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